Simple Combined Model for Nonlinear Excitations in DNA 
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We propose a new simple model for DNA denaturation bases on the pendulum model of 
Englande]>i and the microscopic model of Peyrard et al.^ so called "combined model". The main 
parameters of our model are: the coupling constant k along each strand, the mean stretching y* of 
the hydrogen bonds, the ratio of the damping constant and driven force f/F. We show that both 
the length L of unpaired bases and the velocity v of kinks depend on not only the coupling constant 
k but also the temperature T. Our results are in good agreement with previous works. 

PACS numbers: 87.10.+e, 87.14.-g, 87.15.-v, 87.83. -fa 
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I. INTRODUCTION 

The structures of DNA molecule has been extensively 
studied during the last decade theoretically and experi- 
mentally. DNA molecule consists of a pair of molecules, 
organized as strands running start-to-end and joined by 
hydrogen bonds along their lengths, the diameter of the 
double stranded DNA is 2 nm indeed, but the length 
could be much longer, usually. The separation of double 
helix structure of DNA is an important beginning point 
in the informatics replication process of DNA that needs 
to reproduce the living matter. Accordingly, the problem 
of DNA denaturation has been the subject of extensive 
experimental and theoretical investigation. It is not only 
of practical interest due to its significance in biology but 
also of fundamental interest. 

Englander et al} first constructed a theory of a ge- 
ometrical soliton by carrying the base rotations to the 
point that one side of the helix was given an entire extra 
360° twist. This theory is an explanation of the open 
states of DNA due to the nonlinear excitations along 

On the other hand, a simple lattice model (PB model) 
for the denaturation of the DNA double helix was pro- 
posed by Peyrard et al.,'^'^'^ and more advances are dis- 
cussed in 7.8,9l fl0lfl]| . They introduced the Morse inter- 
action potential that depends on the transverse stretch- 
ing of the hydrogen bonds between the complementary 
base pairs. These models could be used to explain 
the experimental observations on biological phenomena 
of DNA (see Refs. ^M^- 

However, recently it was 
arguedi^ii^ii^iiiii^ that thermodynamic characterization 
of the thermal fluctuations may differ from a dynamical 
characterization, which points to the need for a thorough 
understanding of the dynamical effects in this highly non- 
linear and cooperative material. 

Note that, in Englander model, the strands of a double 
helix are modeled as two chains of pendula (the bases), 
in which one of the strand is dynamic while other strand 
is fixed. By using the continuum approximation, the 
relative motion of the chain has been described by the 



sine-Gordon equation for a single pendulum in the grav- 
itational potential. In Peyrard-Bishop model, only the 
out-off-phase displacements stretch the hydrogen bonds. 
In the continuum limit approximation and thermody- 
namic limit, the out-off-phase motion of the bases is de- 
scribed by an equation which is formally identical to the 
Schrodinger equation for a particle in the Morse poten- 
tial. Therefore the two above models started from the 
slightly similar two-line simple caricatures of DNA (see 
the figures [1] and [2] bellow) describing physical processes 
derivate from the same stretching of hydrogen bonds and 
motions of the similar bases, then both were reformulated 
to onc-dimcnsional motions in some effective potentials 
with slightly similar equations. Because of that it is pos- 
sible to combine the models. 

In this paper we propose a new simple model for DNA 
denaturation based on the combination of two above sim- 
ple models of the DNA, i.e., (E) and (PB) models. By 
using this combined model (EPB model), we study the 
temperature dependence behaviors of the nonlinear exci- 
tations in DNA. 

The paper is organized as follows: Section [ll] gives 
a brief review of Englander model pi A[) and Peyrard- 
Bishop model piBI) . In Sec. IIIIl wc suggest the com- 
bined model by extending their works to consider the 
temperature dependence of the main parameters of the 
model. Some applications are shown in Sec. IIV[ we in- 
vestigated the dependence of unpaired bases length and 
velocity of kinks on the temperature, coupling constant, 
and external driven force. Concluding remarks are given 
in Section IVl 



II. BACKGROUND 

A. Englander pendulum model (E model) 

The model of DNA proposed by Englander et alr^ is 
schematically represented in Fig. [l] where the dynam- 
ics of one of the strands of the DNA is represented as 
a chain of pendula; leaving the other strand fixed (the 
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FIG. 1: Sketch of the Englender models (a) The ground 
state of the double hehx. (b) Kink sohton in the sine-Gordon 
model, the bases of the upper chain twist through 2ti over 
a characteristic length L (the mean number of open bases) 
along the line, (c) The angle 9n is the deviation of the upper 
base of the nth pair with respect to the plane defined by the 
fixed base pairs in the lower strand. 



one at the bottom in Figs. [Iji, [Tja), the base pairs of the 
(upper) strand behave like pendula in an gravitational 
field caused by the tendency of the base pairs of the two 
strands to be facing each other (see Fig. [Dd). It must 
be realized that this model describes only the dynamics 
of the base pairs around the sugar-phosphate backbone, 
which is assumed to be fixed as well. 

Introducing a damping constant 7 and an external 
driven force F, these pendula are described by the 
discrete, dc-driven, damped sine-Gordon model. The 
Hamiltonian of the coupled-pendulum system can be 
written as 



dt 



+ mgJ{l-cosen)+j-^+F9„j, (1) 

where m and I is the common mass and length of pendu- 
lum, 9n and gn is the rotation angle and the gravitational 
acceleration of pendulum at site n, respectively, k is the 
coupling constant along the strand. 

In the mechanical unzipping or denaturation experi- 
ments, one of the DNA strands is attached to a glass 
bead, which is pulled by a glass micro-needle, while the 
other strand is attached to a glass plate, which serves as 
fixed reference point. The glass bead is then pulled at 
constant velocity, opening or unzipping the double chain 
by consecutively breaking the hydrogen bonds of the base 



pairs, and the force used to pull it is recorded. In this 
context, the driven force F in Eq. ([T|) represents the ac- 
tion of the pulling on the end of one strand on the region 
where unzipping is taking place,^^ The damping constant 
7 is introduced to fix the force scale only."^ 

In the case of homopolymer, i.e., gn — g and the pen- 
dulum spacing h in the chain is small, by using the con- 
tinuum approximation 0„(t) — > 9{x,t), we obtained the 
damped sine-Gordon equation 



d?9 



-j*dt9-F*, 



(2) 



here we use the notations 7* — j/mP and F* = F/mP. 
The soliton solutions of the kink type could be obtained 
by solving this equation. 

The physical meaning of a kink solution in the context 
of DNA modeling: the bases of the upper chain perform a 
complete, smooth turn around the sugar-phosphate back- 
bone, from 9n = to 9n = Stt. The kink joins a sector of 
the chain where bases are closed, On — 0, to another one 
where bases have performed a complete turn, 9n — Stt. 
The chain is said to be open in its last part. 



B. Peyrard-Bishop microscopic model of DNA 
(PB model) 

With a particular interest in thermal stability, a mi- 
croscopic model to describe the dynamics of DNA denat- 
uration was introduced by Peyrard and Bishop in 1989j^ 
In Fig. [5] we represent PB model of DNA, where b* and 
a* are the distance between two strands and the spac- 
ing between bases in a strand, respectively. The PB 
model has successfully reproduced the essential features 
of thermally induced denaturation of long DNA chains. 
This model has been used to estimate the melting curves 
of very short heterogeneous DNA segments, in excellent 
quantitative agreement with experimental data.^^-^^ It 
provides the characteristic multistep melting of DNA se- 
quence, and also has been used to investigate the prop- 
erties of dynamic transport in DNA chains. 

For each base pair, there are two degree of freedom 
Un and Vn , which correspond to the displacements of the 
bases from their equilibrium positions along the direction 
of the hydrogen bonds connect the two bases in a pair. 
The potential for hydrogen bonds is approximated by 
a Morse potential V{r) = _D[exp(— ar) — 1]^ with two 



FIG. 2: Sketch of the PB model. 
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parameters D and a. Neglecting the inhomogeneities, 
using a common mass m for bases and the same couphng 
constant k along each strand, the Hamiltonian for PB 
model reads 

HpB = ^^l^miul+vl) + V{Un-Vn) 
n ^ 

+ ^ ^ [(Wn - U„_i)^ + {Vn - Vn-l)'^] | • (3) 

The motions of the two strands can be separated in terms 
of the variables x„ = {u„ + w„)/\/2, ?/„ = (u„ - u„)/\/2, 
which represent the in-phase and out-of-phase motions, 
respectively. Because only the out-of-phase displace- 
ments yn stretch the hydrogen bonds, the component Ti.y 
of Hamiltonian ([3]) with respect to y is taken 



1 



k 



Vn-l 



f + V^(V2y„) . (4) 



As the model is assumed to be homogeneous, the re- 
sult does not depend on the particular site n considered. 
Using the transfer integral (TI) method, with the con- 
tinuum limit approximation and thermodynamic limit 
(large number of bases pairs N ^ oo), the main result is 
dominated by the lowest eigenvalue, so the mean stretch- 
ing {y) of hydrogen bonds is given by 



iv) = {voiv) |y|<Po(y)> = / vliy)ydy 



(5) 



The mean stretching (y) is interesting for the study of 
DNA denaturation, it increases rapidly around a partic- 
ular temperature which is a characteristic of DNA de- 
naturation. The denaturation temperature is not only 
sensitive to the parameters of the hydrogen bonds which 
bind the two strands, but also very sensitive to the in- 
trastrand interaction constant. 



III. COMBINED MODEL (EPB MODEL) 

Comparing the two above models, we see that the Eng- 
lander's model explains the denaturation of DNA as the 
motion of soliton excitations along the chain, while the 
PB model has successfully describe the open states of 
DNA by separating the relative motion of two strands at 
the critical temperature. Therefore, we suggest a com- 
bined EPB model which is based on E-model with includ- 
ing the parameters of PB-model, i.e., taking into account 
the stretching y* of hydrogen bonds 



y 



vliy)ydy, 



(6) 



where the normalized eigenfunction (po{y) is obtained by 
solving the Schrodinger equation for a single particle in 
the Morse potential 



1 



2/32 fc dy 



92 

7rT + ^(2y) 



ip{y) = eip{y), 



(7) 



in which (3 = l/ksT with fc^ being Boltzmann's con- 
stant, then we have 



My) = {V2ay/^ 



(2d)'^-i/2 



X exp 



id- l/2)V2ay 



= exp ( -de 
1)1 ^ 



-\/2ay 



(8) 



Equation ([7]) has a discrete spectrum when d ~ 
{f3/a){kDY^^ > 1/2. In the case T < = 2y/kD/akB 
the states of particle are localized, while they are delo- 
calized for T > Tm. So that the critical temperature 
Tjn is considered as the melting temperature of DNA. 
This equation is also formally identical to the Schrodinger 
equation for a quasi-particle with an effective mass m*: 



(9) 



where mg is the mass of free electron, and Tr = 300 K 
is room temperature. The dimensionless parameter A is 
defined by 



A = R 



H 



{ksTr 



(10) 



in which Rh = 13.6 eV and an = 0.053 nm are the 
Rydberg constant and Bohr radius of hydrogen atoms, 
respectively. Note that A is proportional to the ratio of 
constructive and destructive energies of DNA. 

Note that, the effective mass depends on the temper- 
ature and the structure constant of single strand (see 
Fig. [3|). At room temperature T = 300 K, we found the 
value of effective mass m* ^ 22.87mo. We show in Fig. 3] 
the temperature dependence of the stretching y* of the 
hydrogen bonds for two values of k, with D = 0.33 eV 
and a = 18 nm^-'^. 

Now, we consider the influence of temperature on the 
parameters of the Englander model with the substitu- 
tions: the length and mass of the pendulum I = y* and 
m = m*, respectively, g = g* with g* being the effective 
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FIG. 3: The ratio of m*/mo as a function of temperature T 
with k = 0.2 eV/nm^. 



4 



1.0 

0.8^ 

0.6 

0.4 

0.2 



T{K) 



200 



300 



400 



500 



FIG. 4; The stretching y* of hydrogen bonds as a function 
of temperature T for two values of k: (dashed curve) k — 
0.2 eV/nm^ and (dotted curve) k = 0.23 eV/nm^. 



In the case F = and in the presence of damp- 
ing, 7^0, the only possible value for the velocity is 
V = 0. When both damping and force are present, the 
balance between the two effects leads to kinks propa- 
gating at a constant, nonzero velocity. An analytical 
expression for that velocity can be easilyderived from 
energy-conservation arguments (see Refs.l^fiot). 
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IV. SOME APPLICATIONS 
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FIG. 5: The ratio g* /go as a function of temperature T with 
k = 0.23 eV/nm^, go = 9.81 m/s^, D = 0.33 eV and a = 
18 nm'\ 



gravitational acceleration and defined by 

1 d 
m* ay* 



2Da 



[exp(- 2ay*) - exp{-ay*)] . (11) 



In Fig. [5] we present the ratio g* /go as a function 
of temperature with g^ — 9.81 m/s^. We see that 
g* decrease rapidly to zero at a critical temperature 
Tm ~ 335 K, this temperature is so-called "melting tem- 
perature" and should be discussed in next section. 

It is well known that, in the absence of dissipation and 
force 7 = F = 0, Eq. ^ possesses soliton solutions of 
the kink type, whose expression is 



6±{x,t) = ±4arctan exp 



2h 



'S- t 

v 



v/1 - {v*f 



(12) 



in which, the plus or minus sign stands for kinks or an- 
tikinks, respectively, and < ?; < 1 represents their ve- 
locity, i is a length of unpaired bases such that the mean 
number of open bases along the chain of pendula^ and 
defined by 



2h{k/m*g*y*Y''^. 



We assume that the spacing of bases in a same strand 
does not depend on temperature, because the interaction 
between the bases in each strand is stronger than that 
between two strands. We chose the spacing h — 0.34 nm 
for all numerical calculations. 

For comparing to the results of PB-model^, we used 
the following parameters: D — 0.33 eV, a = 18 nm"-'^. 
By using the numerical calculations, we present in Fig. [6] 
the temperature dependence of the kink velocity for 
three values of coupling constant k ~ 0.2 eV/nm^, 
k — 0.3 eV/nm^, and k = 0.4 eV/nm'^. The velocity 
V = V* /vq increase with increasing the temperature and 
more rapidly close to wq = 1 at the melting temperatures 
T,„ - 323 K (a), ~ 395 K (b), and ~ 458 K 
(c) . These results are in an excellent agreement with the 
results of Ref. i where is about 330 K, 400 K, and 
460 K for three values of k, respectively. The particular 
temperature Tm is a characteristic of DNA denaturation, 
because the hydrogen bonds between the strands are bro- 
ken at this stage. 

Figure [7] shows the temperature dependence of the 
length L for three values of coupling constant k above. 
In the same manner as the kink velocity v, the length 
L increases rapidly at the temperature T ~ Tm, so the 
mean number of the open bases is very large, i.e., the 
base-pair is said to be open at the melting temperature. 
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FIG. 6: Variation of velocity w as a function of temperature 
T with -f/F ^2. (a) k = 0.2 eV/nm^ (b) k = 0.3 eV/nm^ 
and (c) k = 0.4 eV/nm^ 
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FIG. 7; Variation of length L as a function of temperature 
T. (a) k = 0.2 eV/nm^ (b) k = 0.3 eV/nm^ and (c) k = 
0.4 eV/nm^ 




of 0.23 eV/nm^ to obtain a reasonable denaturation 
temperature which are obtained from the experimental 
result»i^ - 330 K to 345 K. Note that our combined 
model works well for long DNA sequences only. It is well- 
known that the finite length effect of DNA sequences is 
important, so that our present results might not appli- 
cable for short chain of DNA, such as in [l2| with 27 
base pairs and 21 base pairs sequences. The influence of 
finiteness length effects of DNA sequences in our model 
will be investigated in the future work. 

Finally, we show in Fig. [10] the dependence of kink 
velocity on the driven force F for three values of tem- 
peratures: (a) T = 340 K, (b) T = 335 K, and (c) 
T ~ 330 K. The velocity increases with increasing the 
driven force and more rapidly at high temperature. The 
melting temperature could be controlled by varying the 
external driven force F. 
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FIG. 8: The dependence of kink velocity on the temperature 
T and coupling constant k with "f/F = 2. 




FIG. 9: The dependence of length L on the temperature T 
and coupling constant k. 

One of important parameters of the model is the cou- 
pling constant, as k increases the denaturation temper- 
ature increases. This is consistent with the increase ob- 
served experimentally in the presence of reagents that 
increase the hydrophobic interactionsi^ Figures [5] and 
|9] show that the values of k must be of the order 



FIG. 10: The dependence of kink velocity on the driven force 
F with 7 = 0.1, k — 0.22 eV/nm^ for three values of temper- 
ature: (a) T = 340 K, (b) T = 335 K, and (c) T = 330 K. 



V. CONCLUSIONS 

In this work we suggest a combined EPB model to de- 
scribe the dynamics of DNA denaturation, which is based 
on Englander's model with parameters derived from mi- 
croscopic PB-model. Using this combined model (EPB 
model) with the motion of bases in both x and y direc- 
tions, we investigated the temperature dependence be- 
haviors of the nonlinear excitations in DNA. We also 
studied the dependence of reasonable denaturation tem- 
perature on the coupling constant of bases and the driv- 
ing force. The temperature dependence of the kinks ve- 
locity and the length of unpaired bases are shown and 
discussed. 

This simple model could be applied to study other 
physical quantities of E-model such as the effects of in- 
homogeneity arising from the genetic sequence and the 
influence of temperature on the kink propagation. These 
are objects of our future investigations. 
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